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Abstract 

Group classification of the generalized complex Ginzburg-Landau equations is pre- 
sented. An approach to group classification of systems of reaction-diffusion equations 
with general diffusion matrix is formulated. 

1 Introduction 

Group classification of differential equations is one of the central problems of 
group analysis. It specifies non-equivalent classes of equations and open the way 
to applications of symmetry tools such as constructing and group generation of 
exact solutions, separation of variables, etc. One of the goals of group classifi- 
cation is a priori description of mathematical models with a desired symmetry 
(e.q., relativistic invariance). 

The first (and very impressive) achievements in group classification belong 
to S. Lie who who had classified second order ordinary differential equations 
and specified all cases when such equations can be integrated in quadratures 
PP. Lie had presented also a group classification of an entire class of partial 
differential equations, namely, linear equations with two independent variables. 
In particular, it was Lie who for the first time describes group properties of the 
linear heat equation 

The next step in classification of heat equations was made by Dorodnitsyn [3] 
who had classified nonlinear diffusion equations 

u t -u xx = f(u) (1) 

where / is a function of u = u(t, x) and subscripts denote derivations w.r.t. the 
corresponding variables 1 . This result was extended by Fushchych , Serov [I] and 
Clarkson and Mansfield jH] to the case of non-classical (conditional) symmetries. 

The results of group classification of equations (0) play an important role in 
constructing of their exact solutions and qualitative analysis of the nonlinear heat 
equation, refer, e.g. to 

In paper [3] a more general equations with nonlinear diffusion were classified which include 
Q as a particular case 



In the present paper we perform the group classification of systems of the 
nonlinear react ion- diffusion equations 



Ut — A(au — v) = f (u, v) 
v t — A(u + av) = f 2 (u, v) 



(2) 



where u and v are function of t, real constant and A is the 

Laplace operator in R m . We shall write (J2J) also in the matrix form: 

U t - AAU = f(U) (3) 

where A is a matrix whose elements are A 11 = A 22 = a, A 12 = —A 21 = b, 

Mathematical models based on equations (J2J are widely used in mathematical 
physics, biology, chemistry, etc. Here we present only two significant examples. 

• The nonlinear Schrodinger (NS) equation in m— dimensional space: 

(id t + A)iP = F(i/>,iJ>*) (4) 

is a particular case of (J2J). If we denote if) = u + iv, F = f\ + if 2 then (jlj) 
reduces to the form Q with A = ( ^ J V 

Equations (j3J) with various nonlinearities F are used in nonlinear optics, non- 
linear quantum mechanics [7] , they serve as one of basic models of inverse 
scattering problem 0. The most popular models are connected with the 
following nonlinearities |9|: 

F = F{if)*if))if), F = {if)*if)) k if), F = (ip*ip)%ip, F = \n(ip*if))if) 
One more interesting particular case of the NS equation corresponds to 

(td t + A)ip = (ip-r) 2 ; 

in this case (JIJ) is a potential equation for the Boussinesq equation for func- 
tion V — d t (if> — if)*). 

Group classification of the NS equation has been performed in paper [TU] . 

• Generalized complex Ginzburg-Landau (CGL) equation 

W T -(1 + if3)AW = F(W, W*) (5) 

also can be treated as a particular case of system (J3J). Indeed, representing W 
and F as W = (u + iv),F — ^(f 1 + if 2 ) and changing independent variable 

t — > t = (3t we transform (jSJ to the form (jHJ) with A = ^ ^ j . The 

standard CGL equation corresponds to the case F = W — (1 + ia)!^^! 2 . 
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Thus the symmetry analysis of equations (J2J) has a large application value 
and can be used, e.g., to construct exact solutions for a very extended class of 
physical and biological systems. The comprehensive group analysis of systems 
((21) is also a nice "internal" problem of the Lie theory which admits exact general 
solution for the case of arbitrary number of independent variables X\, x%, . . . , x m . 

We notice that group classification of equations (J2J) by no means is a standard 
problem of group analysis of partial differential equations which can be solved 
with direct application of well-known algorithms. Because of presence of two 
arbitrary elements, i.e., f 1 and f 2 , this classification needs a rather nontrivial 
generalization of the approach [Hj used for classification of equation (JTJ). 

Equations (JSJ) with arbitrary invertible matrix A were classified in paper jllj . 
To our great a pity, mainly due to typographical errors (made during the editing 
procedure), presentation of results in [TT] was not satisfactory 2 . 

The present paper is the first from the series in which we present the completed 
group classification of coupled reaction-diffusion equations (JHJ) with arbitrary (i.e., 
invertible or singular) matrix A. Moreover, we present a straightforward and eas- 
ily verified procedure of solution of the determining equations which guarantees 
the completeness of the obtained results. 

We also indicate clearly the equivalence relations used in the classification 
procedure, i.e., present explicitly the equivalence groups for for all classified equa- 
tions. In addition, we extend the results obtained in ^1] to the important case 
of non-invertible matrix A and more general equations including both the first 
and second order derivatives with respect to spatial variables. 

Let us note that there are three ad hoc non-equivalent classes of equations 
(JHJ) corresponding to the following forms of matrices A 

where a is an arbitrary parameter. Moreover, any 2x2 matrix A can be reduced 
to one of the forms (JHJ) using linear transformations of dependent variables and 
scaling independent variables in (J3J). 

The NS and CGL equations correspond to matrices A of form /. The general 
equations (JHJ) with such matrices (i.e., generalized CGL equations) are the main 
subject of group classification carried out in the present paper while the cases 
II and /// will be considered in the following publications. Nevertheless till 
an appropriate moment we will consider equations with all types of matrices A 
enumerated in ©. 

2 The tables presenting the results of group classification have been deformed and cut off. 
It is necessary to stress that it was the authors fault, one of whom signed the paper proofs 
without careful reading. 
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2 Determining equations and equivalence trans- 
formations 



In the first stage we restrict ourselves to group classification of equations (JHJ with 
invertible matrix A. Moreover, till an appropriate moment we consider equations 
(J3J) with arbitrary number n of dependent variables. 

Using the standard Lie algorithm (or its specific version proposed in [TT] ) 
one can find determining equations for the functions rj, £ a and ir a which specify 
generator X of the symmetry group admitted by equation (JHJ): 



where a summation from 1 torn and from 1 to 2 is assumed over repeated indices 
v and b respectively, and a temporary notation u = Ui,v = u 2 is used. In a more 
general case of n dependent variables U = column(ui,U2, ■ ■ ■ ,u n ) the repeated 
indices b run over the values 1, 2, • • ■ , n. 

We shall not reproduce the deduction of the determining equations here (refer 
to [TT]) but present them directly. 

Dependence of rj, £ u and n b on U is defined by the following relations: 

Vu a = 0, C Ub = 0, < cUb = 0. (8) 
So from (jHJ i] and £ u are functions of t and and, ir u is linear in u a . Thus: 

7T a = N ab u b + B a (9) 
where N ab , B a are functions of t and x v only. The remaining equations are [TTj : 

2AC U = -5^( Vt A + [A, N]), Vxvt = 0, (10) 



£t - 2AN Xv - AA£" = 0, 

Vtf k + N kb f b + (N t kb - AA ks N sb )u b + B k t - AA kc B c = (B a + N ab u b )f k a . { } 

Here and A are matrices whose elements are N ab and A ab , S ab is the Kronecker 
symbol. 

In accordance with (|9* J — fjllj ) the general form of the related generator (J7J is 

m 

X = \K + cr^G^ + + fiD- (C ab u b + B a )d Ua , , 

+^ x ^d Xv + ud t + Pfi d Xlt 1 ' 

where the Greek letters denote arbitrary constants, B a are functions of t, x, and 
C ab are functions of t satisfying 

C ab A bk _ A ab c bk _ g ( 13 ^ 
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and 

K = 2t(td t + xpdvj - ^{A~ l ) ah u b d Ua - tmu a d Ua , 
D = td t + 

Here A ab and (A~ 1 ) ab are elements of matrix A and matrix inverse to A respec- 
tively. 

In accordance with (jllj) equation (JHJ) admits symmetry operator (|T2|) iff the 
following classifying equations for f 1 and f 2 are satisfied: 

(Xt(m + 4) + n)f* + (fx 2 + + 7 e 7 ^x M ) (A" 1 )^/' 

+C ab / 6 + Cfu b + 5 t a - AA ab B b (15) 

= (5 s + C 3b u b + Amtw, + (fx 2 + + 7 e 7i ^x M ) (A- 1 )**^)/^. 

Thus the group classification of equations (J3J) with a non-singular matrix A re- 
duces to solving equation (JT5j) where A, /x, cr^, u;^, 7 are arbitrary parameters, £> a 
and C afc are functions of (£, x) and t respectively. Moreover, matrix C with ele- 
ments C ab should commute with A. 

We notice that relations (|T2" |) -(fT5 ]) are valid for group classification of sys- 
tems dHJ) of coupled reaction-diffusion equations including arbitrary number n of 
dependent variables U = (ui, 11%, . . . u n ) provided the related n x n matrix A 
be invertible [TT]. In this case indices a,b,s,k in ([T2" ]) -(fT5 ]l run over the values 
1,2... n. 

We will solve classifying equations (fTKj) up to equivalence transformations 
U -> U = G(U,t,x), t -> t = T(U,t,x), x -> x = X(U,t,x) and / -> / = 
F(?7, £, x, /) which keep the general form of equations © but can change functions 
f 1 and / 2 . The group of equivalence transformations for equation (jHJ) can be found 
using the classical Lie approach and treating f 1 and f 2 as additional dependent 
variables. In addition to the obvious symmetry transformations 

£— >t' = t + a, -> x' M = R^ v x v + bp (16) 

where a, fe^ and _R Mi , are arbitrary parameters satisfying R^ U R^\ = S^x, this group 
includes the following transformations 



u a -> + 6 a , / a -> X 2 K ab f b 

t ^ /X . ft ^ /\ x ri 



2+ „ _ w„ (17) 



where i^ afe are elements of an invertible constant matrix K commuting with A, 
A 7^ and b a are arbitrary constants. 

For the case when n = 2 and matrix A belongs to type I © the transformation 
matrix has the following form 

K= (k 2 *? + *^°- (18) 
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It is possible to show that there is no more extended equivalence relations valid 
for arbitrary nonlinearities f 1 and f 2 . However, if functions f l ,f 2 are specified, 
the invariance group can be more extended. In addition to transformations (|17j) 
it includes symmetry transformations generated by infinitesimal operator ()12j) . 
and can include additional equivalence transformations (AET). We will specify 
AET in the following. 

3 Basic, main and extended symmetries 

Thus to describe Lie symmetries of equation © (whose generators have the gen- 
eral form <mj)) ^ is necessary to find all non-linearities f 1 and f 2 which satisfy 
the corresponding classifying equations (fTo^) . To solve these rather complicated 
equations we use the main algebraic property of the related symmetries, i.e., the 
fact that they should form a Lie algebra. In other words, instead of going throw 
all non-equivalent possibilities arising via separation of variables in the classifying 
equations we first specify all non-equivalent realizations of the invariance alge- 
bra for our equations. Then using the one-to-one correspondence between these 
algebras and classifying equations (fT5|) we easily solve the group classification 
problems for equations (JBJ). 

Equation (fT5|) does not include parameters vl/^, v and p v present in (j!2)l thus 
for any f 1 and f 2 equation (jSJ) admits symmetries generated by the following 
operators 

P = d t , P x = d Xxl = x^d Xv - x v d Xfi . (19) 

Infinitesimal operators (j!9|) generate the evident symmetry transformations 
(llfijl which form the kernel of invariance groups of equation 0. For some classes 
of nonlinearities f 1 and f 2 the invariance algebra of equation (JHJ) is more extended 
but includes as a subalgebra. We will refer to (fTT)|) as to basic symmetries. 

Let us specify one more subclass of symmetries of equation (j3J) which we call 
main symmetries. The related generator X has the form (|T2*j) with = v = 
Pv = o v = u u = 0, i.e., 

X = jj,D + C ah u b d Ua + B a d Ua . (20) 

The classifying equation for symmetries (j2*U|) can be obtained from (j!5|) by 
setting /i = o a = uj a = 0. As a result we obtain 

(fi5 ab + C ab )f b + Cfu b + B a t - AA ab B b = {C nb u b + B n )f a Un . (21) 

It is easily verified that operators (|2T)|) and (j!9)l form a Lie algebra which is 
a subalgebra of symmetries for equation (JHl) (this algebra can be either finite of 
infinite dimensional). 
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On the other hand, if equation Q admits a more general symmetry (|12|) with 
o a 7^ or (and) A ^ 0, uj^ ^ then it has to admit symmetry ()20|) also. To 
prove this we calculate multiple commutators of (fT2j) with the basic symmetries 
(119 J) and use the fact that such commutators have to belong to symmetries of 
equation (JHl), i.e., generate their own classifying equation (fT5|) . 

Let equation (|2J) admits symmetry (fT2^) with ^ 0, = p M = v = X = 
uj k = 0, i.e., 

X = a u G l/ + fiD + (C ab u b + B a )d Ub . (22) 
Commuting Y with P„ we obtain one more symmetry 



Y, = -J±(A- l y b u b d Ua + B«d Ua + pPp (23) 



The last term belongs to the basic symmetry algebra (fTTH) and so can be 
omitted. The remaining terms are of the type ([20)1 . 

Thus supposing the extended symmetry (|2*2*j) is admissible we conclude that 
equation (jSJ) has to admit the main symmetry also. 

Commuting (|2*B*|) with Pq and Pa we come to the following symmetries: 

Y, v = B; u d Ua , Y, t = B; t d Ua . (24) 

Any symmetry (|22jl -(|24 jl generates this own system (fl3|) of classifying equa- 
tions. After strait forward but rather cumbersome calculations we conclude that 
all these systems are compatible provided the following condition is satisfied 

(A~ 1 ) ab f b = (A- l T b u h f a Un . (25) 

If (J2SJ) is satisfied equation Q admits symmetry (J2~2~j) with [i = C ab = B a = 0, 
i.e., Galilei generators G v of (JT3J). 

Analogously, supposing that equation (jSJ) admits extended symmetry (|T2*j) 
with A 7^ and u fl = we conclude that it has to admit symmetry (j22j) with 
/i ^ and a v ^ also. The related functions f 1 and f 2 should satisfy relations 
(j25j) and ()21|) . Moreover, analyzing possible dependence of C ab and B a in the 
corresponding relations (fTB]) on t we conclude that they should be ether scalars 
or linear in t, i.e., C ab = fi ab t + v ab . Moreover, up to equivalence transformations 
(fT7|) we can choose B a = and reduce the related equation system (fl~5j) . (j2~Tf to 
the following equations: 

( m + 4 )/a + ^abfb = {fi kb Uf) + mM ,)|g, 

u ab fb + ^ab Ub = u kb Ub dIl W 

where constants v ab and /z ab are non-trivial in the case of the diagonal diffusion 
matrix only. 
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Finally for general symmetry (|T2*|) it is not difficult to show that the condition 
u a 7^ leads to the following equation for f a 



(A- 1 ) kb (f b + 1 u b ) = (A- 1 r b u b ft a . 



(27) 



We notice that relations (|23j) and (}2Tj) are particular cases of (|2*T|) for fi = 0, 
C ab = (A-iyb and ^ = o, C ab = e^iA' 1 )^ respectively Thus if relation fl23) is 
valid then, in addition to (fTljl equation admits the symmetry 



Thus it is reasonable first to classify equations (jHJ) which admit main symme- 
tries ()2())1 and then specify all cases when these symmetries can be extended. 

The conditions when system Q admits extended symmetries are given by 
relations (|23j) - (l2Tj) . 

Results of Section 3 are valid for equations Q with arbitrary invertible diffu- 
sion matrix. In the following we restrict ourselves to the case when n = 2 and ma- 
trix A has the form / given in (JHJ), i.e., to the case of generalized CLG equations. 
In other words, we will classify the following systems of coupled reaction-diffusion 
equations: 



which are particular cases of general systems (j3J) when matrix A has the form / 



4 Algebras of main symmetries 

In accordance with the plane outlined above we start with investigation of main 
symmetries (|2T)j) admitted by equation (jBTIj) . 

The first step of our analysis is to describe non-equivalent Lie algebras A 
of operators ()20)1 which can be admitted by this equation. We shall consider 
consequently one-, two-,..., n-dimensional algebras A. 

For any type of matrix A enumerated in @ we specify all non-equivalent 
"tails" of operators (J20J), i.e., the terms 



These terms can either be a constituent part of a more general symmetry (|2T)|) 
or represent a particular case of ()20|) corresponding to fi = 0. Thus the problem 




u t — A(au — v) = f (u, v) 
v t — A(at> + u) — f 2 (u, v) 



(30) 



N = C ab u b d Ua +B a d Ua . 



(31) 
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of classification of algebras A includes a subproblem of classification of algebras 
of operators (f3*T|). 

Let equation (JHJ) admits a one-dimensional invariance algebra whose basis 
element have the form (f3~T|) . and does not admit a more extended algebra of the 
main symmetries. Then commutators of iV with the basic symmetries Pq and 
P a should be equal to a linear combination of N and operators (fTTJj) . It is easily 
verified that there are three possibilities: 

1. C ab = fi ab , B a = fi a , 

2. C ab = e A '/i afc , B a = e A '/i a , (32) 

3. C ab = 0, B a = e xt+w - x ii a 

where /j, ,fi a , A, and to are constants, and matrix with elements fi ab should com- 
mute with A. In the case when A is of the form (??) constants [i ab are restricted 
by the following relations: /i 11 = /i 22 , /i 12 = — /i 21 . 

To specify all non-equivalent operators (}3~T]) . (J3~2"j) we use the isomorphism of 
(13 1|) with 3x3 matrices of the following form 

00 \ / \ 

B 1 C 11 C 12 ~ fi 1 /i 11 ^ 12 . (33) 

£2 £21 £13 y ^2 _^12 ^11 J 

Equations (J3H|) admit equivalence transformations (fTTj) which change the term 
AT (|3*T|) and can be used to simplify it. The corresponding transformation for 
matrix (|3^j) can be represented as 

g - 0' = L^LT 1 (34) 

where {/ is a 3 x 3 matrix of the following special form 

/ I h l K ] K 2 j . (3- r )) 




Up to equivalence transformations (jBIJ), (f33|) there exist three matrices g, 
namely 



9i = 1 , g 2 = 1 , </ 3 = a -1 . (36) 






In accordance with (|31jh (}3*2*j) the related symmetry operator can be repre- 
sented in one of the following forms 

Xx = fxD - 2(g a ) bc u c d Ub , X 2 = e xt (g a ) bc u c d Ub (37) 
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or 




X 3 = e xt+ "- x d U2 . (38) 

Here (g a )bc are elements of a chosen matrix (}36j) . b, c =0, 1,2, u = column (1, u, v). 

Formulae (|H7|) and (J3%|) give the principal description of one-dimension alge- 
bras A for equation (J3U|) . 

To describe two-dimension algebras ^4 we classify matrices g (}3*3*j) forming two- 
dimension Lie algebras. Choosing one of the basis elements in the forms given in 
(f3Tlj) and the other element in the general form ()33j) we find that up to equivalence 
transformations (f3*4*j) there exist three two-dimension algebras of matrices g 

M,\ = {91,93}, A 2 , 2 = {#2, #4}, ^2,3 = {01,02} (39) 
where gi,g%-,gz are matrices (JHHJ), and 
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Algebras A 2j i and A 2j2 are Abelian while the basis elements of A 2 ^ satisfy 

[91,92] = 92- 

Using (|3*9*j) we easily find pairs of operators ()20|) forming two-dimension Lie 
algebras. Denoting 

e a = (e a )abUbd Ua , a = 1,2 

we represent them as follows: 

< fiD + ei + ute 2 , e 2 >, < fJ-D + e 2 + vte\, e x >, , , 

< fiD - e!,isD - e 2 >, < iqei + G x e 2 , F 2 e x + G 2 e 2 > ^ ' 

for e\,e 2 belonging to algebras A 2) \, A 2t2 , and 

< fiD - ei,e 2 >, < [iD + ei + vte 2 ,e 2 > (41) 

for e x and e 2 belonging to A 23 . 

Here {Fi, Gi} and {F 2 , G 2 } are fundamental solutions of the following system 



F t = \F + vG, G t = aF + 7 G (42) 

with arbitrary parameters A, v, a, 7. 

The list (I40|) - (|4*T]) does not includes algebras of the types 

< Fe, Ge > and < fiD + e vt ^ x e, e vt ^ x e > (43) 
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(with F, G satisfying (021)) which are either incompatible with classifying equa- 
tions ([15)1 or reduce to one- dimension algebras. All the other two-dimension 
algebras A can be reduced to one the form given in (jjpj) , (JiTf using equivalence 
transformations (|17jl . 

In analogous way we can find three- and four-dimensional algebras of operators 
(f2~Uj) . Thus we have two three-dimension algebras of matrices (j33J) 

^3,1 : ei = e 2 = # 2) e 3 = # 4 ; ^ 
^3,2 : e\ = 92, e 2 = g 3 , e 3 = # 4 

and the only four- dimension algebra: 

M,i : ei = s-i, e 2 = flf 3 , e 3 = g 4 , e 4 = gr 2 . (45) 
Algebras (jUj) can be generalized to the following algebras *4 

< fiD- 2ei, e 2 , e 3 >, < D + 2ei + 2vte 2 , e 2 , e 3 >, ^ , 

< D + 2ei + 2z/te 3 , e 3 , e 2 >, < ei, F x e 2 + G x e 3 , F 2 e 2 + G 2 e 3 > 1 ' 

while (}4*4*j) generate the following algebras: 

< /mD - 2e x , e 2 , e 3 >, < ei, Z? + 2e 2 + 2/ite 3 , e 3 > . (47) 
Finally, v4 4 2 generates the following algebras A 

< fiD - 2ei, z/D - 2e 2 , e 3 , e 4 >, < e u e 2 , e» t+J/ - x e 3 , e" t+J/ - x e 4 > . (48) 

The list (j4*U|) - (j4*H|) does not include algebras which have subalgebras (}4*3|) 
(which are incompatible with classifying equations (|21|l). 



5 Classification results 

Using results presented in previous section we easily perform group classification 
of equations ([50]). To do it we solve the classifying equations (|2*T|) with their 
known coefficients C ab and 5 a which are defined comparing (j2*U|) with the found 
realizations of algebras A. 

If equation (|HU|) admits one-dimensional algebra A, i.e., one of algebras (pTTj). 
(I38J1 the related functions J 1 and / 2 have to satisfy the corresponding determining 
equation (J2*T1 ) which define J 1 and / 2 up to arbitrary functions. In the case of two 
dimension algebras whose realizations are given by relations ([40)1. (|4~Tj) we have a 
system of two determining equations corresponding to two basis elements which 
usually define f 1 and f 2 up to arbitrary parameters. In addition, we control 
the cases when equation (|HUJ) admit extending symmetries G M , and K (JHJ), 
i.e., when the found functions f 1 and / 2 satisfy conditions (|2"H|l . (|2*7|) and (|2*Hjl 
respectively. 
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The next important step is to find additional equivalence transformations 
admitted by equations (j3J) with specified non-linearities f 1 and f 2 . These trans- 
formations are relatively easy calculated using the standard Lie algorithm and 
treating f 1 and f 2 as additional variables. 

We will not reproduce here the related routine calculations but present their 
results in Tables 1-3, where non-linearities f 1 , f 2 and the related symmetries are 
specified. Additional equivalence transformations are given here also. 

The list of possible AETs is present in the following formulae 

1. u — > exp(ut)u, v — > exp(ut)v, 

2. u — ► u cos ut — v sin ut, v — > v cos ut + u sin out, 

3. u — ► exp(tut)u, v — > v + cu^, 

4. u — ► u + out, v — > v, 

5. u — > u, v — > v + out, 

6. u — ► exp(uout) (w cos(<ru;£) + vsm(aout)) , (49) 
?; — > exp(z/tut) (ucos(cru;£) — « sin(cro;t)) , 

7. ii — > exp(2u;t) (w cos(cxu;t 2 ) — u sin(cra;t 2 )) , 
f — > exp(2a>i(:) (u cos(cru;t 2 ) + M sin(<7u;t 2 )) , 

8. it — > exp(Aa;t 2 ) (w cos(2cji(;) + t>sin(2a;t)) , 
f — > exp(Acjt 2 ) (u cos(2c<jt) — w sin(2co>£)) 

where the Greek letters denote parameters whose values will be specified in the 
tables. In contrast with (|16|h (JT7j) these transformations are valid only for equa- 
tions with some special non-linearities f 1 , f 2 specified in the Tables. 

Table 1 presents non-equivalent equations which admit the main or the main 
and extended symmetries. The conditions for non-linearities which extend the 
symmetries are specified in the fourth column. The additional equivalence trans- 
formations are presented in the last (fifth) column. 

In Table 2 the non-linearities are collected which correspond to the main 
symmetries only. The additional equivalence transformations are specified in the 
fifth column. 

In Table 3 symmetries of a subclass of equations (f3T)j) are specified. The 
additional equivalence transformations are given in square brackets and placed 
in the fourth column. 

In Tables 1-3 G^, G M and K are operators (fT4*|) were A is matrix of type I 
(jBJ), ^juOe) is an arbitrary solutions of the Laplace equation: 

The Greek letters denote arbitrary parameters, moreover, up to equivalence trans- 
formations we restrict ourselves to e — ±1 and k = 0, ±1. 

In Table 1 R = [u 2 + v 2 ) 1 , z = tan -1 (-) , -Fi and F 2 are arbitrary functions 
of Rexp(fiz). 
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Table 1. Non-linearities and extendible symmetries for equations (J30|) 



No 


Nonlinear terms 


A/Tain 
IVlcLlIl 

symmetries 


!\ r~\ iH l t l An n 1 

/\QQlL10Ilcil 

symmetries 


AFT 

. \ r J i 

(HHD 


I. 


J — ur i -|- C-T2 
— (/iU + v) , 


e Kt (/^ - 8 Z ) 


(j a , if 


2, if /i 




/ 2 = - uF 2 




(j a , if 


= K = 




+KZ (u — UV) 




/i = a, k = 


2. 


p = e uz R a (\u- fMv), 
p = e uz R a \\v + fiu) 


aD — ud u — vd v , 
uD — ud v + vd u 


Gq, if v = aa 

m. 
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In Table 2 and F2 are arbitrary functions whose arguments are specified 
in Column 3, R and z have the same meaning as in Table 1. 



Table 2. Non-linearities and non-extendible symmetries for equations 

(HDD 



No 


Nonlinear terms 


Argu- 
ments 

nf F 
01 r a 


Symmetries 


AET 


I. 


p = u u+1 F l 
p = u u+l F 2 


u 

V 


vD — ud u — vd v 


1 if 

u = 


2. 


p = Pu + F 1 , 
p = - KU + F 2 


V 


e^+^m K (x)d u 


4 if 

k = = 


3. 


p = e KV F u 
p = e KV F 2 


U 


kD — d v 




4. 


p = uiFx + elnu), 
p = v (F 2 + e\nu) 


V 

u 


e £t (ud u + vd v ) 




5. 


f 1 — e Xz (uF\ + vF 2 ), 
p = e Xz {vF 1 -uF 2 ) 


Re uz 


XD + vRd R - d z 


6, a = 1 
if A = 


6. 


p = KV »+\ 
p = /3v u+1 




vD — ud u — vd v , 
*o(x)d v 


4 


7. 


P = ne\ 
P = (3e* 




D - d v , V (x)d u 


4 


8. 


P = /ilnv, 
P = e\nv 




D + ud u + vd v 
+ ((// — ea)t 


4 


9. 


f 1 = \ 
P = In u 




D + ud u + vd v + td v , 
^o(x)d v 


3,5 
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Table 3. Symmetries of equations (J30|) with non-linearities 

f 1 = (fiu — av) hiR + z(Xu — vv), f 2 = (fxv + au) lni? + z(\v + uu) 



No 


Conditions 
for coefficients 


Main symmetries 
and AET (gOJ) 


Additional 
symmetries 


1 


A = 0, 

fi = V 


e^<9 z , e^\Rd R + atd z ) 
[AET 7 if /x = 0] 


G Q if aa = 0, 
/x^0 


G Q if 

a = i/ = 0,ff^0 


2 


A = 0, 


[AET 6 if fj,u = 0] 


G Q if/t^0 
aa = v — fi, 
or a = 0, /x 7^ 


G Q if 

(7,(T = 77 // = D 


3 


fi + v = 2u 
A ^ 


X 3 = e^ (2XRd R + (u-fjL)d z ), 

C z "T (/-A. 3 

[AET 6 if/i + i/ = 0], 
[AET 8 & 1 if /x = i/ = 0] 


G a if cuo 7^ 

Ci^/X — IS ) — Z/\ 


G a if 

az/ = —A, uo = 


4 


A ^ 0, 5 = 1, 

U/-J- — U/Q ZL 1 


e?+ t (\Rd R +(u+-n)d g ), 

[AET 6 if /xi/ = Act, 
& 1 if /x = cr = 0] 


G Q if /x/v 7^ Act, 
A = a(/x — z/ + acr) 


G a if 77/x = Act, 
A = a/x 


5 


5 = -1 


exp(co>ot) [2Acosti?9^ 

+ ((i/ — /x) cos t — 2 sin £) 9 a ] , 

exp(uj t) [2XsmtRdR 

+ ((i/ — /x) sin £ + 2 cos t) 9«] 


none 



6 Discussion 

We perform group classification of generalized CLG equations (130)1 . i.e., find all 
non-equivalent equations of the considered type and describe their symmetries. 
The obtained results can be used to construct exact solutions for equations which 
admit sufficiently extended symmetries, using the standard Lie algorithm |12j . 
The other application is to search for models with a priory required symmetry, 
e.g., Galilei-invariance. 

In Tables 1-3 all non-equivalent generalized CLG equations are listed together 
with their symmetries and additional equivalence transformations. More exactly, 
we specify here only extensions of the basic symmetries (JT§j) and did not consider 
linear equations. 

First we notice that the usual CGL equation appears as a particular case 
of the classification procedure. The related non-linearities and symmetries are 
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present in Table 1, Item 2 when p = 0. In addition to the basic symmetries (|19p 
this equation admits the dilatation symmetry and symmetry ud v — vd u , which 
correspond to scaling of dependent and independent variables and multiplying 
solutions of the CLG equation by a phase factor. 

In accordance with our classification there exist eight types of generalized 
CLG equations defined up to arbitrary functions Fi and F 2 depending on vari- 
ables indicated in the third column of Table 1, Item 1 and Table 2, Items 1-5. 
Nonlinearities corresponding to the most extended symmetries are collected in 
Table 1. In particular there are nonlinearities corresponding to Galilei- invariant 
equations (|3*U|) (refer to Table 1, Item 1 for a = p, v = 0): 

u t — A(au — v) = uFi + vF 2 , 
v t — A(av + u) — vFi — uF 2 

where Fi and F 2 are arbitrary functions of Re az . This system can be rewritten 
as a single equation for a complex function W = u + iv: 

W t -(a + i)AW = TW (50) 

where T — Fx — iF 2 is a complex function of real variable £ = In \W\ + az with 
z being a phase of W. 

The standard CLG equation does not belong to the class (jHOjl and so is not 
Galilei invariant. On the other hand, setting in (|50|) a = we come to the 
Galilei-invariant subclass of the NS equations (jlj. 

The non-linearities enumerated in Table 3, Item 2 of Table 1 and Items 6-9 
of Table 2 are defined up to arbitrary parameters. The most extended symmetry 
is indicated in Item 2 of Table 1 and corresponds to the following equation for 
complex function W = u + iv: 

W t - (a + i)AW = a(e az \W\) p W (51) 

where a = X + ia is a complex parameter and p — — . 

In accordance with the above equation (joTjl admits Lie algebra of the Shrodinger 
group including operators P M , J^ u (fTUj) and also generators of dilatation D, Galilean 
transformations and conformal transformations K (JHJ). Setting in (J5l|) a = 
we reduce (j51j) to the very popular NS equation with critical power 4/m non- 
linearity. 

If p — then equation ()51|) admit all the above mentioned symmetries except 
the generator K of conformal transformations. 

In general we indicate six classes of equations (J30|) which admit symmetries 
Gfj, and so are invariant with respect to Galilei group. Namely, in addition to 
(IHOJ) . (JBTj) there are the following Galilei-invariant equations: 

iW t + AW = -aW\n\W\ (52) 
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which corresponds to Item 1, v = of Table 3, and 

W t - (a + i)AW = cW(\n\W\+az) (53) 

where c is a complex number equal to ia, a{i — a) or fi + ia for versions indicated 
in Items 2, 3 or 4 of Table 3 correspondingly. 

Non-linearities collected in Table 2 correspond to equations (I30j) which admit 
the main and basic symmetries only. 

Thus we present completed group classification of systems of reaction-diffusion 
equations with square diffusion matrix of type / ©• 

The additional aim of this paper is to present an effective approach for solving 
classifying equations ([To]) . It was demonstrated in Sections 3 and 4 that the 
problem of group classification of equations (jHJ) can be effectively reduced to 
searching for the main symmetries (jHJ). We also make a priori specification of 
these symmetries using the fact that they should form a basis of a Lie algebra. 
The idea of such a specification was proposed in papers [T3| and [T3j . 

In our following publications we use the approach described here to classify 
react ion- diffusion equations with diagonal and triangular diffusion matrix. In 
other words we plane to complete classification of equations (JHJ) with general 
diffusion matrix whose non-equivalent versions are given by formulae ©• 
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